function iiriirTest()
  % Load the reference data
  data = load("../data/gse2.txt");
  data = double(data);
  % Set an IIR filter
  b = [0.000401587491686, ...
       0.0, ...
      -0.001606349966746, ...
       0.0, ...
       0.002409524950119, ...
       0.0, ...
      -0.001606349966746, ...
       0.0, ...
       0.000401587491686]';
  a = [1.000000000000000, ...
      -7.185226122700763, ...
       22.615376628798678, ...
      -40.733465892344896, ...
       45.926605646620146, ...
      -33.196326377161412, ...
       15.023103545324197, ...
      -3.891997997268024, ...
       0.441930568732716]'; 
  % Filter data
  len = length(data);
  %data(1:20) = 0.0;
  %data(len-20:len) = 0.0;
  yfilt = filtfilt(b, a, data);
  %H = dfilt.df2(b, a);
  %bn = H.Numerator;
  %an = H.Denominator;
  %yfilt = filter(bn, an, data);
  %yfilt = filter(bn, an, flipud(yfilt));
  %yfilt = flipud(yfilt);
  %mean(data)
  %yfilt = H.filter(data);  yw = flipud(yfilt);
  %yfilt = H.filter(yw); yfilt = flipud(yfilt);
  %yfilt = filter(b, a, data);
  % Write the data
  fileID = fopen('iiriirReference.txt','w');
  for i=1:length(yfilt)
     fprintf(fileID, '%.13e\n', yfilt(i));
  end
  fclose(fileID);
end
